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We present a theory of orbital magnetic dynamics for a chiral p-wave superconductor with broken 
time-reversal symmetry. In contrast to the common Landau-Lifshitz theory for spin ferromagnets, 
the case of orbital magnetism cannot be described in terms of local magnetization density. Hence it is 
impossible to define unambiguously the spontaneous magnetic moment: the latter would depend on 
conditions of its experimental investigation. As an example of this we consider orbital magnetization 
waves and the domain structure energy. 
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Superconductors with unconventional pairing mecha- 
nisms have attracted a lot of attention during the past 
decade. A prominent feature characteristic of uncon- 
ventional superconductivity is the possibility of states 
with broken time-reversal symmetry (TRS), which is ex- 
pressed in the presence of magnetic structure in these ma- 
terials. Broken TRS has been detected in several super- 
conductors, including Sr2Ru04 ZrZn2 ,|2| and UGe2, 
which brought renewed interest to the old-standing 
problem of coexistence of superconductivity and mag- 
netism. 

The layered perovskite Sr2Ru04 is one of the best 
known examples of unconventional superconductors 
with broken TRS. Experimental observations of a 
temperature-independent Knight shift for H _L z (Ref. la) 
and an increased muon spin relaxation below T c (Ref. [l|) 
support the theoretically proposed spin-triplet p-w&ve or- 
der parameter 00,01 A a( a(k) = i d l {k)(a l a v ) a p, with 
d(k) = Aoz(fc 2 , + iky). Such an order parameter has a 
non-zero orbital moment 



l=-i«(k)|kx V k |d i (k))/(d;|d J > ) 



(1) 



which should lead, in principle, to spontaneous magne- 
tization like in usual ferromagnets. However, coexisting 
superconductivity screens out this magnetization, mak- 
ing its experimental detection rather difficult. As a pos- 
sible way to overcome this difficulty, it was proposed 01 
to perform microwave response measurements, where ex- 
citation of spin waves would provide an immediate sig- 
nature for the presence of ma gnet ic order. 

Previous works 0, 0, HH Il2t IhH ] which considered 
magnetic dynamics in superconductors with coexisting 
magnetism (SCFM), assumed a phenomenological model 
in which the magnetization was independent from the 
superconducting order parameter, and its dynamics was 
described by the Landau-Lifshitz equation (LL dynam- 
ics). 14] This model, however, is not so obvious for ma- 
terials such as Sr2RuC>4, where TRS is broken at the 
superconducting transition, and the magnetic properties 
are expected to stem from the orbital part of the multi- 
component superconducting order parameter. Therefore 
a proper treatment of magnetic dynamics in such mate- 



rials should derive it from the dynamics of the supercon- 
ducting order parameter. An analogous derivation has 
been done 15] for the A phase of superfluid 3 He ( 3 Hc- 
A); then a modification of the theory taking into account 
the charge of Cooper pairs would yield magnetic orbital 
dynamics for an isotropic p-wave superconducting elec- 
tron liquid. However, such a derivation cannot be applied 
directly to the case of a superconducting metal because 
of the crystal-field anisotropy. 

In this paper, we address this problem and derive an 
effective orbital magnetic dynamics for a p-wave super- 
conductor in a strong crystal-field potential. The result- 
ing dynamics is not equivalent to the phenomenological 
LL SCFM magnetization dynamics, except for some sim- 
ple cases. Moreover, in general, it cannot be described 
in terms of local magnetization at all! Instead, we show 
that it can be described with the help of a unit vector 1, 
corresponding to the angular momentum of Cooper pairs. 
A similar situation happens in 3 He A phase, where the 
presence of an anomalous term in the current does not 
allow us to express the dynamical equations in terms of 
local angular momentum density. 0, 

We start by considering a Ginzburg-Landau (GL) 
free energy functional for a p-wave superconductor in a 
strong crystal field possessing a tetragonal symmetry D^h 
(which is the symmetry of Sr2Ru04). For simplicity, we 
will assume that the spin does not participate in the dy- 
namics. This assumption is justified, for example, in the 
case when the direction of d(k) is fixed along the crystal 
axis z by strong spin-orbit coupling. The order parameter 
can be decomposed into five odd-parity irreducible repre- 
sentations of D^h symmetry. However, limiting the order 
parameter to a p-wave form and requiring d || z leaves 
only two of them: a two-dimensional T§ — {ik x ,ik y }, 
and a one-dimensional = {zk z }. Having in mind 
the order parameter *~ k x + ik y for Sr2Ru04, we assume 
that the representation is most favorable, having the 
highest transition temperature. It will be described by 
the GL order parameter if = (rj x , rj y ). However, this rep- 
resentation alone is not enough to describe the magnetic 
properties, since it does not allow any magnetic moment 
perpendicular to z, and hence we will need an admixture 



2 



of the higher-energy representation T 1 . The free energy 
then has the form 



I = h + h + fi-s + b 2 /8n, 



(2) 



where b = V x A is the magnetic field. The form of the 
GL functional corresponding to representations T§ and 
is well-known: |18( 



f 5h om=Pi(T)\ff\ 2 + 0i\rf\ 4 + 2 (v* x % - V,V* y ) 2 



+fc\Vx\ 2 \Vy\ 2 



(3) 



for the homogeneous part of the functional corresponding 
to the representation, 

hgrad = K 1 \D x ^ y ■ rf\ 2 + K 2 [\D x rj y \ 2 + |ZVfe| 2 ] 
+K 3 [{D x r, x )*{D v r lv ) + c.c] + K 4 [(D x riy)* D yVx + c.c] 
+K 5 [\D zVx \ 2 + \D zVy \ 2 } (4) 

[where Di = di — i(e* / ch)Ai] for the gradient part of the 
functional, and 



h=P 2 (T)\iP\ 2 + K 6 (\D X ^\ 2 + \D y iP\ 2 ) 
+K 7 \D Z ^\ 2 



(5) 



for the functional corresponding to the T 1 representa- 
tion. For Pi(T) < and P 2 (T) > 0, the ground state is 
determined by fhom5, which gives if = 77(1, ±i), provided 
that 02 > 0, 4/? 2 > 03, and 4 (ft - ft) + 03 > 0. This 
reproduces the order parameter for Sr2Ru04, and it is 
also analogous to the Anderson-Brinkman-Morel (ABM) 
state in superfluid 3 He with the 1 vector parallel to the z 
axis 0| . Assuming the couplings of fshom to be domi- 
nant (consistent with the strong anisotropy and the ten- 
dency to TRS breaking), the order parameter 77 is frozen 
in the ground state 



77 = 77(1, i) and 



V 2 o 



'Pi(T) 



4(ft-ft) + ft' 



(6) 



while the low-lying excitations are described by soft vari- 
ables: ip and the total phase (p. 

The mixing terms between T§ and T± are determined 
by requirements of invariance under the point group, as 
well as under the gauge transformation and the time re- 
versal. The lowest-order homogeneous mixing terms are 



*\2 



(7) 



The gradient mixing terms are taken to second order in 
if). Then there are 2 types of such terms: terms such as 
~ (Dr])*Dtjj and such as ~ (r]Dip) 2 . Writing only those 
terms which do not vanish for the state given by Eq. JBJ , 
we have two terms of the first type: 

[(D z1 p)*B x , y ■ ff+ c.c] © [D* z rf ■ T> Xt y1> + c.c] , (8) 



while for the second type four terms are possible: 
\ff\ 2 \T> x ^\ 2 ffi ri*r)*D x <ipD y ilj + c.c. © 

[ril-rilr [(D x i>) 2 - (D y iP) 2 } + c.c. 



(9) 



Collecting all terms and using Eq. JBJ|, we obtain for the 
free energy: 

f = (Ki + K 2 )\B X ^\ 2 + (K 3 - Ki)i[{D x ri)*D v r, - c.c] 
+2X 5 |^r?| 2 + Qi[{D z iP)*{D x + %D v )n + c.c] + 
Q 2 [(D z r])*(D x - iD v )ip + c.c] + Q 3 i[{D x iP)*D y i> - c.c] 
+Q i {(r ] *) 2 [(D x ^) 2 - (A» 2 ] +c.c.} + 
Q 5 i[(ri*) 2 D x if)Dyip - c.c] + f x + P 3V 2 \^\ 2 + b 2 /8ir. (10) 

As we have already mentioned, the state if) = 0, ff = 
?7(l,i) corresponds to the ABM state with 1 || z. We 
would like to exploit this analogy and express the order 
parameter in terms of the unit vector 1. Since the triple 
(rj x ,i]y,ip) transforms like a vector under rotations, we 
can produce a state with nonzero if) by making a small 
rotation around an axis in the x-y plane. Indeed, 



#«a,n»?CM,0) = 77(1, i, -foe 40 ), 



(11) 



where 6a is a small rotation angle, and <f> the angle be- 
tween the axis of rotation n and y axis. From this 
we identify 7^/77 = — 8a exp[i0] . On the other hand, 
l z 5a exp[i<f)] — (l x + Hy), so we can make a replacement 
ijj = —rj(l x + il y ). Thus the presence of a small admix- 
ture of the order parameter ip is analogous to tilting the 
vector 1 away from the z direction. We stress that this is 
only correct as long as ip is small. Large deflections of 1 
require a breaking of the state given by Eq. © , and then 
the configuration space (77, ip) might include also states 
in which 1 is not defined at all (TRS-conserving states). 
This is a consequence of strong crystal-field anisotropy 
(for a weak anisotropy, like in 3 He, this would not be a 
problem, since there states with defined 1 - namely, the 
ABM states - are well separated by the energy from other 
states) . 

In addition to the tilting of the vector 1, excitation 
of the orbital mode results also in superconducting flow 
with the velocity v = Vip/m* — (e*/m*c)A. We would 
like to mention in passing that since only small deviations 
of 1 from i are considered, no problem of nonzero curl 
of the phase gradient (Mermin-Ho relation ^|) arises 
here. Expressing Eq. I|1U|) in terms of the new variables 
1_l = {lx,ly), and v (in the harmonic approximation), we 
obtain for the free energy 



1 



1 



! = -^Pxy\Vxy\ +^Pz\v z \ + B X | V • lj_ | 

+S 2 |z •Vxl 1 | 2 + 5 3 |c4lJ 2 + C v • V x U 



-C an -v z ■ V x lj_ + Zd x l x d y ly 



IM 2 



b 2 /8n, (12) 
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where all the coefficients may be expressed via those of 
Eq. JTJ3: for example, a = 2[P 2 (T)r] 2 + P 3 -q% C = 
QiVo m *i an d C an — {Qi+Q^r^m* . Without anisotropy 
terms oc a and oc Z, this functional is analogous to the 
energy of the 3 He-A phase, linearized with respect 
to small deviations from the z axis. The free energy de- 
termines the (charge) current density: 



._e*8f_ 

J „ n Jtr T Jm j 

nrf ov 



where 



Jtr 



— {Pxy-Vxy + PzV z 



is the transport current, and 



— (CV x l_L-C an V xy x l x ) 
m* 



(13) 



(14) 



(15) 



is the magnetization current [here V xy = (d x , d„ , 0)]. 

Comparing with the LL SCFM model, |3, llOj one can 
see that it differs from the latter not only by anisotropy in 
the London penetration depth ~ 1/pi and the in-plane 
term ~ Z, but, most importantly, by the presence of 
an anomalous term ~ C an in the magnetization current. 
This term renders it impossible to have a consistent defi- 
nition of magnetization for the system (despite the pres- 
ence of magnetization currents!), which makes a crucial 
difference between our p-wave superconductor and LL 
SCFM's. Indeed, this current cannot be expressed as a 
curl of magnetization vector cV x m, and therefore the 
energy of interaction (m*/e*)v ■ j m cannot be cast in the 
Zeeman form — b m. We emphasize here that the anoma- 
lous term appears not due to the anisotropy, but, rather, 
due to the orbital nature of magnetism in our system. 

The magnetic dynamics is generally described by Am- 
pere's law 



4-7T 

(jtr + jm) = V X b , 



by the London equation 



(16) 



(17) 



where £ = (e*/m*c), and by an equation of motion for 1. 
This motion is described by the equation 



at 



-gz x 



Sf_ 
5h 



-gz x 



df_ 
dl± 



Of 
ddih 



(18) 



which is a precession equation for a unit axial vector, 
with a generalized torque given by the expression in the 
brackets. This equation can be obtained from thermody- 
namic conservation requirements by a general procedure 
used for the derivation of hydrodynamic equations. 
Its generalization for the orbital dynamics of 3 He-A has 
been considered in Ref. [Eq. (|18|) is obtained by as- 
suming small deviations from the equilibrium and zero 



normal-fluid velocity]. Although this equation formally 
looks the same as the LL equation for spin ferromagnets, 
there is an essential difference: the dynamical constant g 
here is an unknown phenomenological parameter, while 
in the LL equation it is given by g — j/Mq, where 7 
is the gyromagnetic ratio, and Mq is the ma gnet ization. 
The value of g was disputed for 3 He-A, [lg, [ljj but its 
determination is a prerogative of the microscopic theory. 
Explicitly the equation of motion for the 1 vector reads 

i_L = -gz x [ol± - 2BiV(V • lj_) - 2B 2 z x V(V x l ± ) z 
~2B 3 d 2 l ± - C£b - C an V x v 2 - zSjdMl - <%)] • (19) 

Equations (|T3*|l - ljT7l) and (|T§|) constitute a full system of 
equations describing the magnetic dynamics. As was ex- 
plained above, the magnetic order parameter 1 cannot be 
identified as a local magnetization, despite being analo- 
gous to it, and hence the dynamics is, in general, more 
complicated than the LL SCFM dynamics. In special 
cases, however, namely in situations where all variations 
are either parallel or perpendicular to the z axis, the 
dynamics can be described in terms of local magnetiza- 
tion, so an effective LL SCFM description is valid (with 
a generalization of using tensorial stiffness parameters 
and the London penetration depth, as well as in-plane 
anisotropy). Still, the value of the effective magnetiza- 
tion determined this way would not be unique, but rather 
different in each case. For example, for magnetization 
waves, propagating in the z direction (perpendicular ge- 
ometry), the free energy is 



f = ^ L \v xy \ 2 -C^b-l ± +B 3 \d z l 
leading for plane waves oc e jq ' 
2B 3 q 2 



a 



'8tt' 



CO 

f) 



± a 



4ir- 



%u}t to a dispersion 

c 2 £V 



(20) 



(21) 



This has a form of a LL SCFM spectrum that was con- 
sidered in Ref. [10| with an equilibrium magnetization 
Mq = £C. Hence the results for the microwave response 
of a LL SCFM given in that paper are directly applicable 
to the case of an unconventional superconductor consid- 
ered here, provided that an effective equilibrium magne- 
tization is chosen as Mo = £C- On the other hand, for 
waves propagating in the x — y plane, the spectrum is 



U3 



( Zq 2 sin^ 



(2B 2 - 



Z sin 2 2<t>^ 2 
o )l 



(C 



Zsivi2cj)^ 2 



47T 



a, 



2 



q 2 + 47rp z £ 2 



(22) 



where (f> is the angle the wave vector q makes with the 
x axis. This expression corresponds to the LL SCFM 
spectrum in the parallel geometry |ll| with a different 
value of equilibrium magnetization, Mq = £(C — C an ) (in 



that analysis no x — y plane anisotropy was introduced, 
so the 4> dependence was irrelevant there). We stress 
again, however, that these results cannot be interpreted 
in terms of a tensorial M . For a general propagation 
direction the dynamics cannot be interpreted in terms of 
local magnetization at all. 

Another issue in which important differences arise 
between our unconventional superconductor and a LL 
SCFM, is the issue of domain walls. A basic feature 
of ferromagnets, both insulating and superconducting, 
which determines the field distribution inside domains, 
is the discontinuity of the magnetic induction Ab z across 
the wall being given by 87rMo, with Mq the equilibrium 
magnetization inside the domains. Moreover, this discon- 
tinuity may be used as an experimental definition of do- 
main magnetization. One might ask, whether in our case 
Ab z /8ir has any meaning of magnetization and whether 
it is related in any sense to dynamic response properties 
discussed above. The answer to the second question is, 
in general, negative. Indeed, in order to determine Ab z , 
one has to examine what happens inside the domain wall, 
which involves a strong deviation from the ground state, 
Eq. ©. For that purpose, the free energy Eq. l(T2*)l is in- 
applicable, and one has to go back to the general expres- 
sion, Eq. J2J. With the energy f^hom being dominant, 
a domain-wall solution may be restricted to the {rj x ,r]y) 
part of the order parameter. For specific parameter val- 
ues, domain-wall solutions have been obtained, [22I Esl ] 
with the discontinuity being determined by the param- 
eters of fsgrad- For example, when /3i,/?3 ^> /?2, and 
K3 = K4, the discontinuity is |22| 



domain 1 d. wall domain 2 



Ab z 



47rPi(T)e*if 3 ^i-^2 



K x + K 2 



(23) 



Thus, Ab z has nothing to do either with the parameters 
C and Can, which determine the magnetization current 
for small deviations, or the dynamical parameter g. Only 
in the limit of very weak crystal anisotropy, when the 
vector 1 is well defined even inside the domain wall as 
shown in Fig. ^ the discontinuity Ab z is determined by 
a combination of C and C an . Indeed, in this case the 
magnetization current inside a domain wall is given by 
the isotropic form 



[CV x 1 - C on l(l • V x 1)] , 



(24) 



and then the field discontinuity for a Bloch wall (where 1 
rotates in the wall plane) is given by Ab z = 87r£(C— C an ), 
corresponding to an effective magnetization £(C — C an ). 
On the other hand, for a Neel wall (where 1 rotates in 
the plane perpendicular to the wall) the discontinuity 
becomes Ab z = 8tt£C, corresponding to an effective mag- 
netization £C. 

As to the first question, Ab z /8ir does have the corre- 
spondence to magnetization in the usual ferromagnets. 



a) 
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FIG. 1: Domain wall between two domains with opposite di- 
rections of the 1 vector in the case of weak anisotropy: a)Bloch 
wall; b)Neel wall. Note that when the anisotropy is strong 1 
is not defined inside the wall. 



To see it, let us recall that for usual ferromagnets the 
domain-structure energy is given by the sum of a local 
domain-wall energy (which may be considered an effec- 
tive parameter independent of the fields) and an electro- 
magnetic energy inside the domains (including the energy 
of Meissner currents and the Zeeman term —b z Mo). In 
our case, the Zeeman term is absent, and instead the 
domain-wall energy contains a contribution due to the 
domain- wall currents. However, for thin domain walls 
this contribution can be transformed to an effective Zee- 
man form: 



jyAy 



dx 



Mb z (x)dx, 



(25) 



where M = Ab z /8ir, and the domains are in the x direc- 
tion. Hence Ab z /8-K may be interpreted as an effective 
magnetization of each domain, which allows the electro- 
magnetic part of the domain-wall energy to be cast in 
the usual Zeeman form, so that the remaining part is 
entirely local (i.e., independent of the field distribution). 
This static "magnetization" determines the field in the 
domains, but is unrelated, in general, to the dynami- 
cal response of the material. As an application of this 
result, one can immediately obtain a criterion of stabil- 



5 



ity against the formation of domains for our material. 
For this purpose we can use the result of Krey's analysis 
[24| for a LL SCFM, that the uniformly magnetized con- 
figuration becomes unstable when W < 4irAM§, where 
W is the local domain-wall energy (calculated excluding 
electromagnetic effects), and A the London penetration 
depth. The same result is valid for our case if M defined 
above is substituted for Mq, and the in-plane penetration 
depth is used, A~ 2 — * 4ir£, 2 p xy . Of course, W should be 
calculated separately, as was done in Refs. Eil andE^l 

In summary, we have studied the orbital magnetic dy- 
namics in a p-wave superconductor with strong crystal- 
field anisotropy. The dynamics is essentially different 
from the Landau-Lifshitz dynamics for a superconducting 
ferromagnet. The most important difference is that the 
directional order parameter 1 (orbital moment of Cooper 
pairs) does not lead to a definite spontaneous magnetiza- 
tion (magnetic moment per unit volume). While in sim- 
ple cases one can introduce an effective magnetic-moment 
density similar to that in the Landau-Lifshitz dynamics, 
the value of this density varies from case to case. As ex- 
amples of these cases we have considered magnetization 
waves along and normal to the main crystal axis, and the 
energy of the domain structure. 

This work was supported by the Israel Academy of 
Sciences and Humanities. 
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